A phase-space approach to directional switching in semiconductor ring lasers 
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We show that a topological investigation of the phase space of a Semiconductor Ring Laser can 
be used to devise switching schemes which are alternative to optical pulse injection of counter- 
propagating light. To provide physical insight in these switching mechanisms, a full bifurcation 
analysis and an investigation of the topology is performed on a two-dimensional asymptotic model. 
Numerical simulations confirm the topological predictions. 
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I. INTRODUCTION 

Semiconductor ring lasers (SRLs) have received in- 
creasing attention in recent years owing to their pos- 
sible applications in photonic integrated circuits. This 
has lead to both theoretical and experimental investiga- 
tions, ranging from fundamental studies of its nonlinear 
dynamical behavior to multiple practical applications 1 - 
ITT] . In particular, SRLs attract a lot of interest due to 
the bistable character of their directional mode operation 
and the possibility to encode digital information in the 
direction of emission. 

Recently, all optical flips-flops based on two coupled 
microrings have been fabricated which could be switched 
by injection of a counter-propagating signal [5]. On the 
other hand, a device based on a single ring would benefit 
of a smaller footprint and could take advantage of the 
recently demonstrated switching schemes based on one- 
side injection |10j . However, single ring devices present 
a reduced stability, which leads to the appearance of un- 
desired dynamical regimes such as alternate oscillations 
[TJ, or noise-induced mode hopping [TT] . 

A full understanding of the dynamical stability prop- 
erties of a SRL as well as an in depth knowledge of the 
phase space structure would allow for robust devices and 
smarter switching schemes which rely on the topological 
properties of the system [10 . A stability diagram of SRLs 
has been proposed, which is based on an asymptotic two- 
dimensional model [12 , including analytical expressions 
for the local bifurcation curves. However, the important 
issue of global bifurcations remained unaddressed. 

The aim of this paper is twofold. First, we provide 
a complete bifurcation analysis of the two-dimensional 
SRL model, which includes all the global bifurcations and 
the codimension-two points that act as organizing cen- 
ters for the dynamical behavior. The phase space topol- 
ogy of the system is systematically investigated for time 
scales which are slower than the relaxation oscillations 
[TZ] . In particular, we identify the role of the different 
physical parameters on the global shape of the invari- 
ant manifolds. Secondly, we propose two novel topology- 
based schemes to switch the operation of the ring from 
clockwise (CW) to counter-clockwise (CCW) mode. Such 
schemes are based on externally steering the system in 



the phase space and are alternative to the standard ap- 
proach of injecting an optical signal counter-propagating 
to the lasing mode [HI [14]. More specific, we will show 
that it is possible to switch in a robust way by current 
and phase modulation. 

The paper is structured as follows. In Sect.|H] we start 
by reviewing the existing models that has been shown 
to adequately describe the dynamical behavior of SRLs 
[T[|12) and we revisit previous results on injection induced 
switches in SRLs [TO] [TSJ [14] . In Section III we provide 
insight in the full bifurcation picture of a two-dimensional 
SRL model, including the global bifurcation lines which 
were not addressed before. In Section [TV] we study the 
shape of the invariant manifolds of the saddle points in 
reduced SRL system and we discuss the effect of the in- 
dividual parameters on the topology of the system. The 
particular shape of these invariant manifolds allows us 
to devise novel switching schemes through current mod- 
ulation and phase modulation, which we discuss in Sect. 
|V) Our predictions are verified numerically based on the 
full rate equation model. Finally, in Sect. |VI| we draw 
conclusions of our analysis and discuss possible future 
work. 



II. SEMICONDUCTOR RING LASER MODELS 

We consider a SRL operating in a single-longitudinal, 
single-transverse mode. In order to describe the dynami- 
cal behavior of this SRL, we will focus on a rate-equation 
approach allowing for an easy and fast numerical imple- 
mentation. Problems involving wavelength changes, how- 
ever, fall outside the scope of such rate-equation models, 
and a traveling wave model would be more suitable to 
tackle such issues [T51 ITS] . 

In the limit of small outcoupling from the ring cavity, the 
total electric field oscillating in the ring can be written 
as the sum of two counter-propagating waves: 



E(z, t) = Ei{t) exp [i (uj t — k z)] 
+E 2 (t) exp [i (Lu t + k Q z)] + c.c. 



(1) 



Here fco is the longitudinal wavenumber and ojq is the 
optical frequency of the mode. In the slowly varying en- 
velope approximation, the amplitudes E\2 vary on time 
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scales which are slower than 2-k/ujq. The rate-equation 
model is formulated mathematically in terms of two rate 
equations for the slowly varying amplitudes i? 12 and one 
rate equation for the carrier number TV. The equations 
read [D1GJ]: 

E lft = k(1 + ia) [N (1 - s|^ 2 | 2 - c\E 2A \ 2 ) - l] E 1>2 
-ke^Et,!, (2) 
N = 1 [»-N-N(l-s\E 1 \ 2 -c\E 2 \- 2 )\E 1 \* 

-N(l~ S \E 2 \ 2 -c\E^)\E 2 \% (3) 

where dot represents differentiation with respect to time 
t, k is the field decay rate, 7 is the carrier decay rate, a 
is the linewidth enhancement factor and \x is the renor- 
malized injection current with /1 » at transparency and 
fi sa 1 at lasing threshold. The two counter-propagating 
modes are considered to saturate both their own and 
each other gain due to e.g. spectral hole burning effects. 
Self- and cross-saturation effects are added phcnomeno- 
logically and are modeled by s and c. For a realistic 
device cross saturation is stronger than self-saturation. 

Contrary to the case of solid-state lasers [17], in the 
devices under study, the standing-wave pattern has a 
spatial period much smaller than the carrier diffusion 
length. Therefore, longitudinal variations of the carrier 
density are washed out by the diffusion. As a result, 
the dynamics of such a carrier grating can be neglected 
[5J[THl[n5]. However, a linear coupling between the two 
counter-propagating modes remains. Reflection of the 
counter-propagating modes occurs at the point where 
light is coupled out of the ring cavity into a coupling 
waveguide and can also occur at the end facets of the 
coupling waveguide. These localized reflections result in 
a linear coupling between the two fields characterized by 
an amplitude k and a phase shift (f>k [2DJ . 

In order to gain a deeper insight in the switching mech- 
anism, we have proposed an asymptotic simplification 
of Eqs. (|2])-(|3]) [12 j. A similar reduction has been pro- 
posed to describe the polarization switching behavior in 
VCSELs [5TJ. On timescales slower than the relaxation 
oscillations, it can be shown that the total intensity is 
conserved: 

\E 1 \ 2 + \E 2 \ 2 = n-l >0. (4) 

The slow timescale dynamics is then described by the 
time evolution of two auxiliary angular variables: 

9' = —2 sin <pk sin ip + 2 cos <fcf- cos ip sin 9 

+ J sin 9 cos 9, (5) 
cosdi/)' = aJ sin 9 cos 9 + 2 cos sin ?/j 

+ 2 sin 4>k cos t/j sin 9. (6) 

where 9 = 2 arctan v / j£ 2 l7piF - tt/2 e [-ir/2,n/2] 
represents the relative modal intensity and ijj € [0,27r] 
is the phase difference between the counter-propagating 
modes. Prime now denotes derivation to the slow 
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FIG. 1. Stability diagram for a fixed value a = 3.5 showing 
the different stable operation regimes for J > 0. The different 
black solid lines indicate the different bifurcation currents, 
and are depicted in a polar plot (J, <^). The dotted lines 
depict four cuts for four different values of (j>k, for which we 
show the corresponding Pi-curves in Fig. [2] The eight solid 
grey dots denote the different operating points of which we 
plot the invariant manifolds in Figs. [3] |4j [5] Finally, the 
denotations 1-6 correspond to (1) the Pitchfork bifurcation of 
the OPSS (IPSS), (2) Hopf bifurcation of the OPSS (IPSS), 
(3) Saddle-Node bifurcation the OPSS (IPSS), (4) Sub-critical 
Hopf bifurcation the OPAS (IPAS), (5) Glueing bifurcation, 
(6) a Fold of Cycles, and (7) a double homoclinic bifurcation. 
Multiple Takens-Bogdanov (TB) co-dimension two points are 
also depicted. 

timescale r = kt. Finally, in this reduced model the 
pump current has been rescaled as 

J= k(c- s)(/i- l)/k. (7) 

As the phase space of Eqs. (|5|-([6| is restricted to two 
dimensions, it allows for a clear physical picture of the 
influence of all parameters on the dynamical evolution of 
the variables in a plane. Remark that further simplifi- 
cations to a one dimensional model, as was possible in 
some cases for VCSELs [22j [23] , are not valid for SRLs 

mm- 

Previously [TD], we have reported on the possibility 
of directional mode switching in semiconductor ring 
lasers through optical injection co-propagating with 
the lasing mode. This mechanism is counter intuitive 
as the mode receiving the least energy from the op- 
tical injection prevails after the injection pulse, but 
can be explained due to the folded shape of the stable 
manifold of a saddle point S that is present in the system. 

The understanding of this novel feature in ring lasers 
is thus largely based on the specific structure of the two- 
dimensional asymptotic phase space, in particular the 
shape of the invariant manifolds of the saddle point. In 
this manuscript, we wish to extend this principle by in- 
troducing new ways to switch the SRL between both its 
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directional modes. Exploiting specific phase-space struc- 
tures, we will demonstrate the possibility of directional 
switching through current and phase modulation. In or- 
der to be able to adequately design these schemes, it is 
necessary to have insight into how the phase space pic- 
ture changes with the relevant parameters. Therefore, in 
what follows, we will first perform a full bifurcation anal- 
ysis of the reduced model, which afterwards will provide 
us with the optimal parameter regions for an in-depth 
study of these different switching schemes. 



III. BIFURCATION ANALYSIS 

One advantage of the asymptotic reduction presented 
in Ref. [T^j is the possibility to have an analytical expres- 
sion for all the steady state solutions of the system. The 
second advantage is the two-dimensionality of its phase 
space, which enables not only an analytical study of all 
local bifurcations, but also simplifies the numerical con- 
tinuation techniques necessary for obtaining global bifur- 
cations. Furthermore, it allows for an easy visualization 
of the dynamical behavior in phase-space. 

Two different types of stationary solutions exist. On 
the one hand, solutions corresponding to bi-directional 
emission (6 = 0) exist for all model parameters and are 
identified by their value of tp. For tp — and tp — tt, we 
can distinguish between the in-phase symmetric-solution 
(IPSS) and the out-of-phase symmetric solution (OPSS). 
On the other hand, asymmetric solutions correspond- 
ing to unidirectional emission can emerge from either a 
pitchfork bifurcation of the symmetric solutions or from 
a saddle-node bifurcation and are given by: 



' a sin (p k + cos <p k \ 
; — — tan tp 



9{tp) = arcsin 

\ a cos (pk — sin (pk 
= 2 esc 0{ip) sec 9(ip) 



J OA) 



[sin (pk sin tp — cos 4>k cos tp sin 9(ip)) 



(8) 
(9) 



If the asymmetric solutions originates from a pitchfork 
from the IPSS (OPSS), we will refer to them as IPAS 
(OPAS). This pitchfork bifurcation current can be easily 
obtained in an analytical form: 



±2 



cos <pk + a sin (pk ' 



(10) 



where the upper (lower) sign again corresponds to the 
IPSS (OPSS). We will study the stability of the system 
in a two parameter plane defined by J and (pk with a fixed 
a. If we choose a polar representation where x = J cos (pk 
and y — J sin (pk , the pitchfork bifurcation denoted by O 
in Figure [l] is a straight line x + ay — ±2. For a certain 
parameter regime, the IPAS (OPAS) are born in a saddle 
node bifurcation or fold bifurcation before the pitchfork 
bifurcation has occurred. This Saddle-Node bifurcation 
is given analytically by: 



cos ip 
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(11) 
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FIG. 2. Bifurcation diagrams of Eqs. |5|-l(6| depicting the ex- 
tremes of 6 vs. injection current J. The steady state values of 
9 are denoted by full lines, while the extrema of periodically 
oscillating 9 are indicated with dashed lines. Black (grey) 
color is used for stable (unstable) fixed points or limit cy- 
cles. The dotted lines represent the different bifurcation cur- 
rents, and are denoted as explained in the text. Black (grey) 
OPSS/IPSS/OPAS denotations correspond to stable (unsta- 
ble) OPSS/IPSS/OPAS. From (a) to (d), we have 4> k = 0.2, 
1.08, 1.5 and 2.3, respectively, a = 3.5. 



and it is denoted by © in Figure [T] 

The symmetric solutions can also change stability 
through a Hopf bifurcation at J = Jui'pk) of the symmet- 
ric solutions (OPSS, IPSS) as described by the expression 
below: 



J H = ±4cos( 



(12) 



where the upper (lower) sign corresponds to the IPSS 
(OPSS). In Fig. [l] these Hopf bifurcations lines are seg- 
ments of the circles [(x±2) 2 +y 2 — 4] denoted by ©. These 
Hopf bifurcations lines end at the pitchfork bifurcations 
in a Takens-Bogdanov point with symmetry (TB2) [2~3] , 
The asymmetric solutions can change stability through 
a subcritical Hopf bifurcation given by 



tan 2 



COS (pkOL 



X 

cos(pk + sin (pka / 
sin 2 (pk — cos 2 (pk — a sin 2(pk 
3 sin 2 



(13) 



bk + cos 2 (pk — a sin 2<pk 

These Hopf bifurcations (© in Fig. [T]) emerge from a 
Takens-Bogdanov point (TBi) at the Saddle-Node bifur- 
cation line. 

As previously mentioned, the dynamics is organized 
by two co-dimension two Takens-Bogdanov points. It is 
known that from these points there should also be global 
bifurcations emerging |25j . Therefore, we have comple- 
mented our analytical analysis of the bifurcations in the 
reduced system, with a numerical approach using the 
continuation software package AUTO [26]. From TBi, 
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the point where the Saddle-Node bifurcation line © and 
the Sub-critical Hopf line © meet , emerges a double ho- 
moclinic bifurcation line ©. From TB 2 , the point where 
the Super-critical Hopf bifurcation line © and the Pitch- 
fork bifurcation © meet, emerges a double homoclinic 
connection with a glueing bifurcation © and a fold of cy- 
cles ©. Finally, also the Super-critical Hopf bifurcation 
and the Fold of Cycles coincide in a co-dimension 2 point 
called a Bautin bifurcation point (not shown) [2T1 |2"5"] . 

As the Takens-Bogdanov points involve local bifurca- 
tions, we have obtained analytical expressions. The crit- 
ical phase <Pk B2 °f this Takens-Bogdanov (TB2) point 
with symmetry can be expressed analytically as follows: 



sin 



b TB 2) = 



1 a 
2 



jTB 2 = ±4cos 



1TB2 



(14) 
(15) 



We point out that Figure [T] provides a complete pic- 
ture of the different dynamical regimes for a — 3.5, but 
only for positive values of J, meaning that c > s, being 
the physically relevant parameter region. The bifurca- 
tion lines for negative values of J, however, are exactly 
the same, but due to the particular invariance below of 
the reduced system (Eqs. 5][6) the stability is reversed: 

<t>k-^ 4>k + 7T, (16) 
-J, (17) 
i-> -t. (18) 

While performing a bifurcation analysis, both ana- 
lytically and numerically, many different dynamical re- 
gions of operations can be found in the reduced system 
Eqs. as can be seen from the different Pi-curves 

in Fig. |2J Three different regions of operation A,B,C can 
be identified in Fig. [2]ja-d) according to the relative mag- 
nitude of the linear and nonlinear coupling terms between 
CW and CCW mode. 

At the threshold current, J « 0, laser action starts. 
When operating close to the threshold, at low power, 
the nonlinear coupling between CW and CCW is negli- 
gible when compared to the linear coupling induced by 
the backscattering. The resulting operation is bidirec- 
tional with CW and CCW operating cither in-phase or 
out-of-phase according to <pf.. (see region A in Figure 
[2| . For large values of the pump parameter J, the non- 
linear gain saturation becomes the dominant coupling 
mechanism between the modes and two out-of-phase (in- 
phase) asymmetric OPAS (IPAS) solutions become pos- 
sible. The optical output power is mainly concentrated 
in one propagation direction, called unidirectional oper- 
ation. Because of the device symmetry two OPAS exist: 
one where \Ei\ 2 > \E 2 \ 2 and vice versa. In this regime, 
the device exhibits bistability (See Regime C in Figure 
[21). For intermediate values of J, the linear and nonlinear 
coupling are comparable (Region B) and the dynamics of 
the system as well as the shape of the P-I curve depend 
on the value of <pk ■ A SRL aimed to application such as 
all-optical information storage is expected to be operated 



at the edge between region B and C. 
Depending on the linear coupling phase <f>k, the dynami- 
cal behavior that the SRL exhibits at the transition be- 
tween bidirectional and unidirectional operation, origi- 
nates from instabilities of either the out-of-phase or in- 
phase solutions. The stable bidirectional solution can be 
destabilized, immediately going to stable unidirectional 
operation [see Fig. |2jb)], but there can be many other 
possibilities of which we show a couple of examples in 
the other panels of Fig. [2j E.g. in Fig. |2ja), there exists 
a region where both the stable bidirectional solution as 
the stable unidirectional solution coexist. In Fig. (2jc), 
the destabilisation of the bidirectional solution leads to 
region B, characterized by stable alternate intensity os- 
cillations between the two counter-propagating modes, 
which has been observed experimentally [T] . At this Hopf 
bifurcation point, a limit cycle representing a dynamic 
competition between the two counter-propagating modes 
appears. Furthermore, a combination of the above situ- 
ations is also possible. For instance, in Fig. [2jd) , there 
exists tristability between the stable IPSS and the stable 
OPAS, with the presence of different unstable limit cy- 
cles separating the different basins of attraction of these 
stable solutions. 



IV. PHASE-SPACE REPRESENTATION 

In this section we investigate the phase space (9, ijj) for 
a SRL operating in the bistable regime. In this bistable 
unidirectional regime, Eqs. ([5])-([6]) have four stationary 
solutions: an unstable in-phase (out-of-phase) bidirec- 
tional state in (0,0) [(0, n)]; two symmetric stable states 
CW (clockwise) and CCW (counter-clockwise) at -0 w n 
(tp 0) both corresponding to undirectional operation, 
and a saddle point S in (0,7r) [(0,0)] which is the unsta- 
ble out-of-phase (in-phase) bidirectional solution. As we 
are interested in operating regimes which are relevant for 
optical information storage, we consider here parameter 
values such that we are at the edge between region B 
and C (see Fig. [2j. In particular, we focus on parameters 
such that we are relatively close to the Fold of Cycles 
bifurcation line ©. 

We are interested in the qualitative shape of the invari- 
ant manifolds of the saddle-point S and their dependence 
on the principal parameters of the system. In particu- 
lar, we focus on the stable manifold of S separating the 
basins of attractions of the CW and CCW mode, shown 
in black and white in Figs. [3] [3] [5] (The parameter sets 
under study in these figures are depicted by grey dots in 
Fig.Q. 



1. Influence of J 

We study first the dependence of the phase space struc- 
tures on the pump parameter J [see Eq. Q]. As de- 
fined in Sec. ITT1 J measures the ratio of nonlinear and 
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FIG. 3. Stable manifold at <f>k = 1.5, and a — 3.5 for 
three increasing different values of the normalized current J: 
J = 0.71, 0.739, 0.8, 0.93, from (a) to (d). The basins of at- 
tractions of the CW and CCW mode are given by the black 
and white region. 

linear coupling between CW and CCW modes. Con- 
sider here a) — (1.5,3.5); when varying the param- 
eter J, one observes three different regions of operation: 
bidirectional operation (A), bidirectional operation with 
alternating oscillations (B) and unidirectional operation 
(C) [see Fig. [2] as it has been experimentally observed [T] . 
The shape of the stable manifold of S, and the basins of 
attraction of CW and CCW are depicted in Fig. [3j The 
basins of attraction of CW and CCW fold into each- 
other as the stable manifold of S spirals inwards. As 
J increases, the stable manifold unfolds and the basins 
of attraction separate. Therefore, a large bias current /j, 
or a small coupling amplitude k correspond to unfolded 
basins of attraction. 



2. Influence of cf>k 

It is very hard to determine the value of <pk for a 
SRL device as this parameter describes the phase of the 
backscattering that arises due to different phenomena: 
localized backscattering at the output coupler region, at 
the end facets of the output waveguides and distributed 
backscattering originating from e.g. side- wall surface 
roughness in the ring itself. Therefore, it is useful to 
gain insight in the effect of the linear coupling phase <f>k 
in the entire region from zero to 27r. In order to study 
the influence of <pk on the shape of the manifold, we will 
focus on the region having the stable OPAS at higher 
current (<pk G [0, 7r]). The region (f>k S [it, 2tt] is exactly 
the same if one interchanges everywhere in- phase (IP) 
and out-of-phase (OP). 

One could wonder whether this folded shape of the stable 
manifold is only present if one has a region of alternating 
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FIG. 4. Stable manifold at a = 3.5 for different val- 
ues of the linear coupling phase <j>k'- from (a) to (f) 
J) = (1.5,0.739), (1.64,0.972), (1.78,1.260), (1.36,0.625), 
(1.22,0.625), (1.08,0.636). The basins of attractions of the 
CW and CCW mode are given by the black and white re- 
gion. 

oscillations in the PI scan, as was the case in Fig. [3] We 
have checked the evolution of the manifold shape in the 
region (j>k £ [0, (j\ B2 \ where one does not have a Hopf 
bifurcation [See Fig. kld)-(f)]. 0^ B2 defines the phase 
4>k where the second Takens-Bogdanov bifurcation point 
lies. 

For <p k g [0, 0fc S2 ], one can notice in Fig. ^d)-(f) that 
although one does not observe oscillations in a PI scan, 
the manifold can still have a folded shape. It unfolds 
continuously with decreasing values of (f>k- 
In the region > (j>^, B2 , many different bifurcation sce- 
narios can take place in a PI scan. However, the bifur- 
cation scheme always ends in a Fold of Cycles where a 
stable limit cycle (originating from either the OPSS or 
the IPSS) touches the unstable limit cycle formed in the 
glueing bifurcation. As becomes clear from Fig. [IJa)-(c), 
the manifold unfolds itself with increasing values of 4>k- 
Due to the fact that the distance between the glueing bi- 
furcation and the fold of cycles also increases, the stable 
manifold is already more unfolded by the time the fold 
of cycles occurs. 
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FIG. 5. Stable manifold at {<j> k , J) = (1.64, 0.972) for different 
values of the linewidth enhancement factor a: (a) 3.3, (b) 
3.7. The basins of attractions of the CW and CCW mode 
are given by the black and white region. 



3. Influence of a 

Also the value of the line-width enhancement factor a 
is difficult to determine for a practical device |27j . How a 
affects the dynamical system [Eqs. (|5|-(|6|] is illustrated 
in Fig. [5] We have plotted the manifold shapes for differ- 
ent values of the linewidth enhancement factor a, while 
keeping (<p k , J)=(l. 64, 0.972). We find that increasing 
values of a lead to an unfolding of the stable manifold 



havior is observed. 



for <f) k > fpl" 2 , while for <p k g [0,</> 



^ B2 ], the reverse be- 



This can be explained with the analytical expression 



for the Pitchfork bifurcation of Eq. (10 1. This Pitchfork 
bifurcation is a straight line in Fig. [1J with a slope given 
by — 1/a. Thus, on the one hand increasing the value of a 



for 4> k > (f> 



TB-2 



will shift all the bifurcation lines further 



away from the operating point ((j) kl J) = (1.64,0.972), 
such that the manifold unfolds itself. On the other hand, 
increasing the value of a for 4> k 6 [0, (j>^ B2 ], the pitchfork 
bifurcation will come closer to the chosen (<f> k , J) point 
in the bistable regime, such the manifold will spiral more 
into itself. 

In conclusion, in this Section we have revealed a non- 
trivial shape of the basin of attraction of the system 
which is not observed in other optical systems such as 
VCSELs [H]. The standard approach to switches, con- 
sisting in injection of counter-propagating signal Q21 03] 
correspond, in phase-space, to steering the system very 
close to the final counter-rotating state. However, since 
Figs. |3][4] suggest that other regions of the basin of attrac- 
tion are easier or more useful to target, we investigate this 
possibility in more depth in the next Section. 



V. ALTERNATIVE SWITCHING 
MECHANISMS 

In this section, we propose different switching mecha- 
nisms which are based on the structure of the phase space 
of the system. Consider for instance a SRL operating in 
the CCW mode. In order to switch operating direction, 
one needs to steer the system to the basin of attraction 
of the CW state with an external perturbation such as a 
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FIG. 6. Numerical simulations of Eqs. (p|-([3|, where the cur- 
rent is modulated around fjibias = 1-65 (Juas = 0.739). The 
current modulation Ap is continuously changed in time by 
0.044m with a period of T = 100ns, where m is an integer de- 
scribing the m th modulation step. The insets 7, //, III depict 
the evolution of the system projected in the two-dimensional 
phase plane (8, ijj) for three values of the current modulation, 
respectively for m — 4, 12 and 19. (j> = 1-5) k = 0.44ns -1 , 
a = 3.5, c = 0.01, s = 0.005, k = 100ns" 1 , 7 = 0.2ns" 1 . 



pulse in the bias current or an external optical field. One 
can engineer switches in the following way. Consider the 
laser operating stationary in unidirectional regime, for 
instance CW (Ocw, ''Pew), an d introduce a perturbation 
in the system which shifts the system away from this CW 
state (Ocw, ^Pcw)- It is clear from the previous discus- 
sion that it is possible to use such perturbation to drive 
the system into the basin of attraction of the opposite 
state. Once the perturbation is removed, the system will 
deterministically relax to CCW operation. The switching 
induced by the injection of a resonant optical signal has 
been investigated in 021 031 03] and we will not discuss 
it further here. We propose to induce a switch by: 

1. modulation of the pump current /i, 

2. external modulation the linear coupling phase cj> k . 



A. Current modulation 

Knowing that the OPAS separate for increasing val- 
ues of the normalized pump current J, one possibility 
to induce a switch is by increasing the pump current for 
a certain time, such that the system relaxes to the new 
OPAS state with higher intensity. If one then decreases 
the current fast enough, the system will find itself in the 
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basin of attraction of one of the original lower intensity 
OPAS modes. It will follow the deterministic trajecto- 
ries of this last current value and relax to the targeted 
OPAS state. Which one depends on the amplitude of 
the current modulation determining where in the stable 
manifold one ends up. 

A confirmation of this behavior can be seen in Fig. |6| 
where the full SRL rate equation model [Eqs. (|2])-(|3f] 
is simulated for fibias = 1-65 (J — 0.739), cj) = 1.5, 
k = 0.0044, a = 3.5, c = 0.01, s = 0.005, k = 100ns" 1 , 
7 = 0.2ns . Furthermore, a time- varying modulation 
A/i(m) is applied with a period T = 100ns. A/i(m) 
is given by 0.044m, where m is an integer describing 
the m th modulation step. This A/x corresponds to 
A J = 0.05m. In Fig. [6] one can clearly see that there 
exist regions of A/j, inducing a switch, and others that 
do not, confirming the relevance of the two-dimensional 
stable manifold. 

The insets /, II, III in Fig. (JgJ depict the evolution of 
the system projected in the two-dimensional phase plane 
(9, tp) for three values of the current modulation, re- 
spectively for m — 4, 12 and 19. These insets demon- 
strate that depending on your initial condition (deter- 
mined by An), you either switch to the opposite mode 
or not. More importantly, one can see that the invari- 
ant manifolds of the saddle point as obtained from the 
reduced 2D model are followed by the system described 
with the full SRL rate equation model, confirming the 
validity of the asymptotically reduced system. However, 
one should note that this good correspondence between 
the full rate-equation system and the reduced one is only 
obtained after a certain transient time, determined by 
the damping rate of the relaxation oscillations (ROs). 
We have indeed already mentioned before that the two- 
dimensional system [Eqs. (|5])-(|6])] is only valid on a slow 
time scale (r = kt) , such that the faster ROs should have 
first damped out before one can adequately predict the 
evolution of the SRL with the reduced model. 
In Fig. [6] we have taken an infinitely fast rise and fall time 
for the current modulation, which is not realistic from an 
experimental point of view. In order to understand the 
effect of a finite rise/fall time, we have investigated the 
case of Afj, — 0.0528 (Inset II in Fig. [6]) for different values 
of the fall time. Our simulations (not shown) reveal that 
fall times which are longer that a critical tf a u ss 1.35ns 
do not result in a switch, and the system returns to the 
original state by following the instantaneous equilibrium 
states. 



B. Phase modulation 
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FIG. 7. Numerical simulations of Eqs. ((2j — (J3j , where the lin- 
ear coupling phase is modulated between 4>k,i = 0.7 and (a,c) 
4>k,f = 1.5, (b,d) 4>kj — 4.3 with a period T = 200ns. (a,b) 
depict a timetrace demonstrating the possibility to switch the 
SRL between the two OPAS. (c,d) show the projection of (a,b) 
in the phase portrait of the system in a solid black line. CW 
and CCW are the stable unidirectional states of the system; 
S is the saddle. The stable (unstable) manifold of S is in- 
dicated with a grey solid (dashed) line for the two values of 
the phase tj>k,i and (j>k,f, respectively in dark grey and light 
grey, fj, = 1.63 (J = 0.716), (j> = 1.5, k = 0.44ns" 1 , a = 3.5, 
c = 0.01, s = 0.005, K = 100ns" 1 , 7 = 0.2ns" 1 . 



An important drawback of using current modulation to 
induce directional switches in a SRL, is that it should be 
feasible to electrically modulate the pump current with 
fall times of the order of ns. Modulating the linear cou- 
pling phase could help to circumvent this issue, since the 



modulation could be done optically. Moreover, it is clear 
from Fig. [TJ that one can induce larger effects by mod- 
ulating the phase as compared with current modulation 
and therefore it is more flexible. In what follows, we 
propose two switching schemes that are both based on 
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steering induced by phase modulation. Consider for in- 
stance injection of a non resonant external optical signal 
in the coupler. Far from the resonance condition, the 
main effect of such signal would be to alter the refrac- 
tive index of the coupler without altering the field inside 
the ring cavity. In first approximation, for a standard 
semiconductor, we can disregard the changes in k and 
consider only changes in fa. 

In Fig. [7] we show numerical simulations of the full 
rate-equation system [Eqs. Q-Q], where the linear cou- 
pling phase is modulated between fa, L = 0.7 and fa t — 
1.5 [Fig.[7ja),(c)l, fa j = 4.3 [Fig.0;b),(d)] with a period 
T = 200ns. Fig. [7^a)-(b) depict a timetrace demonstrat- 
ing the possibility to switch the SRL between the two 
OPAS. Fig. [7jc)-(d) show the projection of Fig. [^-(b) 
in the two-dimensional phase space of the reduced system 
in a solid black line. CW and CCW denote the stable 
unidirectional states of the system; S is the saddle. The 
stable (unstable) manifold of S is indicated with a grey 
solid (dashed) line for the two values of the phase faj and 
faj, respectively in dark grey and light grey. The other 
parameters are kept fixed: u = 1.63 (J = 0.716), cj) = 1.5, 
k = 0.0044, a = 3.5,c = 0.01, s = 0.005, k = 100ns" 1 , 
7 = 0.2ns" 1 . 

One can distinguish between two possible methods 
to switch the laser between two OPAS (IPAS) states 
through phase modulation. The first one is based on 
changes in <pk smaller than tt (Acpk < it) and do not 
largely affect the variable ip. Therefore, they involve 
only the areas of phase space which are close to the ini- 
tial modal phase difference ip = ir (ip = 0). The second 
method takes advantage of a phase modulation Acpk > tt 
which affects the complete two-dimensional phase-plane 
(0,ip). While the first method depends on the folded 
shape of the stable manifold of S as discussed in Section 
IV the second one does not and has a operating region 
which can be ten times as large. 

Method 1 is shown in Fig. [7ja)-(c). Starting from 
the clockwise (CW) OPAS solution for fa ti = 0.7 
(OPAS^i.cw), we change the linear coupling phase to 
4>kj = 1-5, such that OPAS^i^cw now finds itself in 
the basin of attraction of the counter-clockwise (CCW) 
OPAS state (OPAS^f^cw) for the dynamical phase pic- 
ture with 4>kj — 1-5. In this new situation, the system 
follows first the stable and then the unstable manifold of 
the saddle point S^f, such that it finally relaxes to the 
stable OPAS^f^ccw state. If we then remove the mod- 
ulation of the coupling phase, the system feels again the 
manifolds of the original phase-space picture (4>k,i = 0.7). 
Since the OPAS^f^ccw state finds itself in the basin of 
attraction of the CCW OPAS state for fa :i = 0.7, it 
relaxes to this OPAS^^ccw state, such that we have 
accomplished a successful switch. 

Since we have seen in Sect. |IV| that the spiraling shape 
of the stable manifold is quite sensitive to parameter 
variations, the necessity of this folding shape for the di- 
rectional switching to occur can be a drawback of this 
method. Therefore, we propose a second scheme, which 



is less demanding on the particular phase space struc- 
ture of the reduced dynamical system. In Fig. [7jb)-(d), 
we again prepare the system in the OPAS^i cw state 
with (f>k.i — 0.7, but now we change the coupling phase 
to faj = 4.3. In doing this, the OPAS^i^cw state is now 
located in the basin of attraction of the in-phase asym- 
metric solution IPAS<pf t cwi such that it relaxes to that 
state. Furthermore, this IPAS^f.cw state finds itself in 
the basin of attraction of the counter-propagating mode 
OPASs i.ccw- By removing the modulation (fa ti = 0.7), 
the system thus follows first the stable and then the un- 
stable manifold of the saddle S^i to end up having suc- 
cessfully switched to the CCW mode OPAS^^ccw- It 
is important to note that this second method works for 
4>kj G [4.27 — 4.47] (for this particular cp^j and current 
J), which is a more than ten times bigger operating re- 
gion than for method 1. 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

To summarise, in this paper we have presented a com- 
plete stability analysis of a SRL which is based on an 
asymptotic two-dimensional reduction of a rate-equation 
model [12 . In particular, we have found the global bifur- 
cation curves corresponding to a homoclinic bifurcation 
of limit cycles and fold of cycles [See Fig. [I] which are 
responsible for the global shape of the invariant man- 
ifolds of the system. Four different Takens-Bogdanov 
codimension-2 points have been located, around which 
the dynamics of the system is organized. 

For parameters which correspond to directional bista- 
bility, the shape of the basins of attraction has been in- 
vestigated. The basin of attraction of the CW and CCW 
are separated by the stable manifold of the saddle S, and 
fold one into each other. The role of the different physi- 
cal parameters of the system in the folding/ unfolding of 
the manifold has been addressed. 

After our topological analysis, we have proposed 
two novel schemes to switch the operation of the ring 
from clockwise (CW) to counter-clockwise (CCW). Such 
schemes relying on the folding shape of the invariant 
manifolds of the system, are alternative to the stan- 
dard approach of injecting a counter-propagating signal 

rjiE3[n]. 

The first of such schemes consists of a modulation in 
the bias current, which drags the system outside the basin 
of attraction of the initial state and releases it when in- 
side the basin of attraction of the final state. Our simu- 
lations reveal the existence of a limiting fall time t f a u for 
the applicability of the scheme. Although the switching 
speed, obtainable with this scheme, is of course limited 
by the electrical bandwidth of the modulation source, all- 
optical methods relying on phase space steering |28j can 
be applied to increase the performance. 

A second scheme has been proposed which consists 
of modulating the backscattering phase <pk- We discuss 
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switching takes place in two different ways when the am- 
plitude of the phase modulation is larger or smaller than 
7T. Numerical simulations have demonstrated the appli- 
cability of the scheme in both cases. 

As a final remark, we want to state that this phase- 
space topology is common for optical systems, for in- 
stance a nonlinear ring resonator described by the Ikeda 
map |29) . Moreover, our reduced two-dimensional SRL 
model can — in the region close to the bidirectional so- 
lution (6 w 0) — also be seen as an optical prototype 
of the symmetric FitzHugh-Nagumo equation, which de- 
scribes excitable systems, and is important in the study 
of neurons [SU) EH] • More generally, we expect our two- 
dimensional system to be relevant for planar dynamical 
systems which are .^-invariant |25j , such that the tech- 
niques discussed here could be applicable to more and 



diverse dynamical systems. 
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